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RESUMEN 

' Presentamos modelos de evolucion qui'mica para el disco de nuestra galaxia. 

■ Tambien presentamos una nueva dcterminacion de X , Y , y Z para M17, una 
, region H ii de nuestra galaxia rica cn elementos pesados. Comparamos nue- 

stros modelos del disco galactico con las abundancias de las regiones H ii. El 
^ . valor predicho por nuestro modelo para AF/AO es muy similar al valor que 

\ obtenemos por medio de las observaciones de Ml? y la abundancia primordial 

^ ■ de helio, Yp. A partir de M17 y Yp obtenemos que ISY j I\Z = 1.97 ± 0.41 

. resultado que concuerda con dos determinaciones de AY/AZ, obtenidas a 

I partir de observaciones de estrellas enanas K de la vecindad solar, que corre- 

Qh' sponden a 2.1 ± 0.4 y 2.1 ± 0.9 respectivamente. Nuestros modelos ajustan 

Q . razonablemente bien el valor de 0/H con el que se formo el Sol. 

to ; ABSTRACT 

We present chemical evolution models for the Galactic disk. We also present a 
new determination of X, Y, and Z for M17 a Galactic metal-rich H II region. 
We compare our models for the Galactic disk with the Galactic H ii regions 
abundances. The AY/ AO ratio predicted from the galactic chemical evolution 
i model is in very good agreement with the AY/ AO value derived from M17 

' and the primordial helium abundance, 1^, taking into account the presence 

■ of temperature variations in this H ii region. From the M17 observations we 
^ i obtain that AY/AZ = 1.97 ± 0.41, in excellent agreement with two AY/AZ 

' determinations derived from K dwarf stars of the solar vicinity that amount 

■ to 2.1 ± 0.4 and 2.1 ± 0.9 respectively. We also compare our models with 
^ , the solar abundances. The solar and Orion nebula 0/H values are in good 

' agreement with our chemical evolution model. 

■ Key Words: galaxies: abundances — galaxies: evolution — H II regions — 
ISM: abundances — ISM: individual: M17, Orion nebula — Sun: abundances 

1. INTRODUCTION 

The main purpose of this work is to study the evolution of the helium 
abundance with respect to the heavy elements as a function of time and posi- 
tion in the Galacti c disk. For this pur pose we will use the Galactic chemical 



cn 
> 



evolution model bv lCarigi et al.l (j2005l ) that has been successful in explaining: 



the observed O/H and C/H abundance gradients in the interstellar medium. 



^Institute de Astronomia, UNAM, Mexico. 
^Centre for Astrophysics, UCLan, UK. 



1 



2 



CARIGI & PEIMBERT 



ISM, the present gaseous distribution in the Galactic disk, the current star 
formation rate, the stellar mass as a function of the Galactic radius, and most 
of chemical properties of the solar vicinity. 

To have useful observational constraints we need accurate X, Y, and Z 
determinations or at least AY / AZ determinations to compare the models 
with observations. In this paper we recompute the X, Y, and Z values for 
the H II region M17, the best Galactic H ii region for which it is possible 
to compute an accurate enough Y value, for this purpose we make use of 
the best observations available and the new He I atomic data needed for the 
helium abundance determination. We also compare the chemical evolution 
model with the AY/AZ determination derived from K dwa rf stars of the 



solar v icini ty with metallicities simi lar or higher than solar by iJimenez et al 
|2003[ ) and lCasagrande et all |2007h . 

We also compare our model with the initial solar O /H value and with the 
Orion O/H value. The initial solar O/H value is representative of the ISM, 
4.5 Gyr ago when the Sun was formed, and the Orion nebula O/H value is 
representative of the present day ISM. 

The model together with the primordial helium determination is also used 
to provide an equation between the Y enrichment and the O enrichment of 
the ISM. This equation can be used to provide the initial Y values for those 
stars for which we can derive their initial oxygen abundances. These initial Y 
values provide meaningful initial abundances for a set of stellar evolutionary 
models with different heavy elements content. 

We adopt the usual notation X, Y, and Z to represent the hydrogen, 
helium and heavy elements abundances by mass, respectively. Based on our 
models we study the increase of helium, AK, as a function of the increase of 
C, O, Fe, and Z by mass. 

In Section [2] we discuss the general properties of the chemical evolution 
models, we discuss inflow models for the Galaxy with two sets of stellar yields 
and present the chemical abundances for the disk at 5 galactocentric distances. 
For the two models discussed we present the increase of helium by mass AY , 
relative to the increase of carbon, oxygen, iron and heavy elements by mass, 
AC, AO, AFe, and AZ. We also discuss the evolution of the Y and O 
abundances for our models, and present an equation that predicts for the 
Galaxy the Y enrichment as a function of the O enrichment of the ISM. In 
Section [3] we present a new determination of the helium abundance for the 
metal rich H ii region M17; this abimdancc is compared with our Galactic 
chemical evolution models. 

In Section 3] we compare abundances for the Orion nebula with those 
of B stars of the Orion association. In Section [5] we discuss the absolute 
calibration of the O/H ratio in the local ISM, this ratio is one of the most 
important observational constraints for Galactic chemical evolution models, 
the absolute calibration of the O /H ratio is obtained based on the Orion nebula 
and the solar abundances. In Section [5] we compare t he solar initial helium 
abundance, inferred from the standard solar models bv lBahcall et al.l (|2006f ). 
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with our Galactic chemical evolution models considering the time since the 
Sun was formed and the presence of gravitational settling and diffusion. The 
conclusions are presented in Section [71 

Throughout this p aper we will us e the p rimordial helium abundance by 



mass, Yp, derived by iPeimbert et alj ([20071) based on direct helium abmi 



dance determinations of metal poor extragalactic H ii regions that amounts 
to 0.2477 ± 0.0029. This re s ult is in excellent agreement with the Yp deter- 
mination by bunklev et all |2008h that amounts to 0.2484 ± 0.0003. This 
determination is based on the il.i,h^ value derived from WMAP observations, 
the assumption of standard big-bang n ucleosynthesis, and the n eutron life- 



time, Tn,, of 885.7 ib .8 sec obtained bv [Arzumanov et al.l (|2000l ). Following 



Mathews et al.[ (|2005[ ) the value of Yp derived from WMAP is revised down- 
wards to 0.2468 ± 0. 0003 by adopting the t„ = 878.5 ± 0.8 sec derived by 
Serebrov et al.l (|2005[ ). and to 0.2475 ± 0.0006 by adopting for t„ the new 
world a verage that amouii t s to 8 81.9 ±1.6 sec, average tha t includes the re- 
sults by Arzumanov et al. ( 2000l ) and Serebrov et al. ( 2005h . 



2. CHEMICAL EVOLUTION MODELS 

2.1. Model Parameters 

We present chemical evo lution models for the Galactic disk using the 
CHEMO code ([Carigi I [l99j) that considers the lifetime of each star until 
it leaves the main sequence. The models have been built to reproduce the 
present gas mass distribution and the present-day O /H values for each galac- 
tocentric dista nce, the values were o btained from observations of H ii regions 
in the galaxy ([Esteban et al.[ [20051 ) . Specifically, the characteristics of the 
models are: 

i) An inside-outside scenario with primordial infalls but without any type 
of outflows. The infall rate as a function of time and galactocentric distance 
r is given by INFALL{r,t) = A(r)e-*/^'"»° + B(r)e-(*^l'='^'')/^<^-^ where the 
formation timescales are Thaio = 0.5 Gyr and rdisk = 6 -f (r/r© — 1)8 Gyr. We 
assume the location of the solar vicinity is r© = 8 kpc. The constants A{r) 
and B{r) are chosen to match, first, the present-day mass density of the halo 
and disk components in the solar vicinity, 10 and 40 Mq^c"^, respectively, and 
second, to rep roduce the radial p rofile total mass in the Galaxy, Mtot{r) = 



50e-('-8)/3-5 (Fenner et al. 200: 



ii) 13 Gyr as the age of the models, the time elapsed since the beginning 
of the formation of the Galaxy. 



iii) The Initial Mass Function (IMF) proposed bv [Kroupa. Tout fc Gilmore 
([l993[ ). in the mass interval given by 0.01 < m/M0 < M.^p, with Mup = 80 
and 60 Mq. This IMF is a three power-law approximation, given by IMF 
a m"" with a = -1.3 for 0.01 - 0.5 Mq, a = -2.2 for 0.5 - 1.0 Mq, and 
a = -2.7 for 0.5 - Mup. 

Note that, our models were computed assuming an IMF with Mjou, = 0.01 
Mq. Kroupa et al. (1993) truncated their IMF at 0.08 Mgbecause they 
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Fig. 1. Newly formed mass of a given element by massive stars, in Mq , ejected to the 
ISM. The initial hea vy elements of the stars amount to Z — 0.02. Co ntinuous lines: 
high w ind yields by iMaeder (|l992l ), dashed lines: low wind yields bv iHirschi et al, 
( 20051 ). 



considered only stars, but in our work we are assuming a non negligible amount 
of substellar objects (0.01 < to/Mq < 0.08). In models with M^p = 60 Mq, 
the mass of objects with m < O.OSMq is 12 % of the total Mgtars, that 
includes stars and remnants; this percentage is practically independent of 
Mup. Even at present, the fraction of mass in substellar objects is unknown 
and we consider that our predicted percentage might be realistic. 

Due to the uncertainties in the current Algasij"), Alstarsij"), and SFR(r) 
values we cannot discriminate between chemical evolution models assuming 
^low = 0.01 MQand Miow = 0.08 Mq. The first ones predict smaller fractions 
of massive stars and LIMS per single stellar generation. Therefore Mio^ = 
0.01 MQmodels with identical galaxy formation scenario, galactic age, and 
stellar yields to those of Miow = 0.08 Mgmodels, require a higher SFR to 
match the present-day 0/H(r) values. The Miow = 0.01 MQmodel with a 
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more efficient SFR predicts a lower Mgas and similar chemical abundances. 
Such model with higher SFR and lowe r Mgas is also able to reproduce the 
observational constraints (jCarigi Ill996( ). 

iv) A star formation rate that depends on time and galactocentric distance, 
that varies from almost constant and low (at large r values) to bursting and 
high (at short r values), this SFR has been represented by the following re- 
lation SFR{r,t) = iyM^ai{r,t) [Mgas + Mstars)°-'^{r,t), in order to reproduce 
the current /H gradient and the gas mass distribution of the Galactic disk 
(|Carigi 1119961 ) , where i/ is a constant in time and space that is chosen in order 
to reproduce the present-day radial distribution of the gas surface mass den- 
sity. A 1/ value of 0.016 is require d when the h i gh-wi nd yields and M„p = 80 
M0are adopted, the best model of lCarigi et al.l (j2005l ). while v values of 0.015 
and 0.010 are required when the low-wind yields with M^p = 60 and M^p ~ 80 
M0are adopted, respectively. 

v) Two sets of stellar yields. Since the main difference between these 
sets is the assumed mass-loss rate due to stellar winds by massive stars with 
Z = 0.02, we will call them high-wind yields (HWY) and low-wind yields 
(LWY), see Figure 1. 



The HWY set is the one considered in the best model (model 1) of lCarigi et al 



(|2005l ). The HWY set includes: A) For ma ssive stars (MS), t hose with 
8 < m/Mg < 80, the following y ields by: a) Ichiefh fc Limongil |2002h for 
Z = 0. 0: b) iMevnet fc Maedeil (|2002l ) for Z = 10"^ and Z = 0.004; c) 
Maedei] (119921) for Z = 0.0 2 (high mass-loss rate yields presented in his Table 



6); d) IWooslev fc Weaver I (|l995l ) only for the Fe yields (Models B, for 12 to 
30 Mq; Models C, for 35 to 40 Mq; while for m > 40 Mq, we extrapolated 
the m = 40 MqFc yields). B) For low and intermediate mass stars (LIMS) 



those wit h 0.8 < m/Mm < 8, we have used the yields by iMarigo et al 



(199 



3 



1998h andlPortinari et al. Hl998l) from Z = 0.004 to Z = .02. G) For Type la 



SNe we have used the yields bv [Thielemann et al. I (|l993l ). We have assumed 
also that 5 % of the stars with initial masses between 3 and 16 M© are binary 
systems which explode as SNIa. 

In the LWY set we have updated t he yield s of m assiv e stars only for Z ^ 
and Z = 0.02 assuming the yields bv iHirschil ([2007|) and lHirschi et al.l (|2005h 
respectively. The rest of the stellar yields are those included in the high wind 
set. 

The main differences between the LWY set and the HWY set are due to the 
contribution of massive stars at Z = 0.02. Therefore in Figure 1 we compare 
the He, G, and O yields for Z = 0.02. The main difference between the HWY 
and the LWY models is due to the stellar yields assumed for massive stars at 
high Z. The HWY assu me a relatively high mass-loss rate for massive stars 
with Z = 0.02 (yields bv lMaedeiill992n . while the LWY assum e a relatively 



low m ass-loss rate for massive stars with Z = 0.02 (yields by iHirschi et al 



20051 ). These difference between a high and a low mass- loss rate produces 



opposite differences in the C and O yields (see Figure 1), the reasons are the 
following: a) a high mass-loss rate produces a high loss of G and consequently 




Fig. 2. Evolution of some common properties for all models of the Galactic 
disk: gas mass surface density (MqPC~'^), star formation rate and infall rates 
(M^pc'^Gyr^^), at four galactocentric distances, 16, 12, 8, and 4 kpc (continu- 
ous, long-dashed, short-dashed, and dotted lines, respectively). 



a high C yield, and b) since C is needed to produce O, the high loss of C 
reduces the O yield. 

Since the solar vicinity and the Galactic disk contain stars and H il regions 
of a broad range of metallicities, our galaxy is a proper laboratory to study 
the AY/AXi behavior at high Z values and to try to observationally test the 
predictions of the HWY models and the LWY models. 



2.2. Results 



The models presented in this paper reproduce the present stellar and gas 
mass distributions in the Galactic disk, the current star formation rate as 
a function of the Galactic radius, the 0/H gradient evolution inferred from 
PNe, the SN rates, the distribution of G-dwarf stars as a fimction of [Fe/H], 
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Fig. 3. Chemical evolution models for the Galactic disk and the solar vicinity {r — 8 
kpc). The left panel shows the C/O evolution in the ISM of the solar vicinity with 
O/H. The right panels show the present-day ISM abundance ratios as a function of 
galactocentr i c dist ance. Continuous lines: high wind yields with Mup = 80 Meby 
Carigi et al. I l|2005l l , dashed and dotted lines: low wind yields with Mup = 60 and 80 
Mq, respectively (t his paper). Filled circ les: H ll regions, gas plus dust values; the 
gaseous values from Esteban et al. | |2005l') have been cor rected for the dust fraction. 
Fille d squares: Dwarf star s from lAkerman et al. (2004). Open circle: Solar values 
from lAsplund et al. I (|2005l ). 



the infall rate , and the evolut i on of [X j /Fe] y s [Fe/H]. See I Allen et al.l (|l998l ). 
Carigi I (|l994l) . ICarigT1 Jl996l ). ICarigTI (j2000l) . 



In Fig. 2 we present some of those model properties, like gas mass surface 
density, star formation rate, and infall rate for four galactocentric distances: 
16, 12, 8, and 4 kpc. Infall rate follows the insidc-outside scenario, the inner 
parts of the Galaxy arc formed faster than the outer parts. Since gas mass 
comes from infall, mainly, Mgas reflects the inside-outside scenario and the 
SFR shows similar behaviour as the gas mass, due to the SFR is proportional 
to the gas mass. 
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TABLE 1 

PRESENT DAY VALUES FROM THE GALACTIC DISK MODELS 



Galactocentric distance O(t/w)(10-^) AF/AC AF/AO AF/AFe AF/AZ 





High wind yields and M^p = 


80 Mr. 






4 


8.89 


6.38 


5.91 


14.11 


1.85 


8 


6.68 


7.01 


4.01 


18.36 


1.67 


12 


4.26 


8.89 


3.29 


21.60 


1.62 


16 


1.95 


10.05 


3.23 


25.28 


1.65 




Low wind yields and M^p ~ 


60 M0 






4 


12.81 


7.13 


4.06 


14.13 


1.67 


8 


6.78 


7.73 


3.79 


18.36 


1.67 


12 


3.28 


9.44 


3.99 


21.52 


1.78 


16 


1.34 


10.18 


4.39 


26.11 


1.90 



In Figure 3 we show O and C gradients in the Galactic disk and the 
evolution of the C/O-O/H relation in the solar vicinity predicted by three 
models that combine different yields and IMF M^p values. Models that assume 
LWY fail to reproduce the C/0 gradient and the C/0 values for halo stars or 
disk stars. The LWY model with Mup = 60 Moreproduces poorly the C/0 
gradient and the C/0 values in disk stars, and predicts C/0 values for halo 
stars higher than observed. The LWY model with M„p = 80 M0does not 
reproduce at all the C/0 gradient, matches partially the C/0 values of disk 
stars, and explains the observed C/0 values for halo stars. The HWY model 
with M„p = 80 M0 reproduces successfully the C/0 Galactic gradient and the 
C/0 values in the solar vicinity. 

Based on Figure 3 we conclude that the LWY model with Ad^p — 60 
MQreproduces the main behavior of C/0 vs 0/H in the solar vicinity but 
cannot reproduce the C/0 Galactic gradient. On the other hand, the HWY 
model reproduces very well the C/0 vs 0/H relation in the solar vicinity and 
the C/0 Galactic gradient. Since the LWY model with M„p = 80 Moproduces 
the poorest fit to the C/0 and C/H observed values it will not be considered 
further. 



Gibson et al I (12006) cons i dering the O yields bv lArnettI (|l99lh . that are 



lower than those by iMaeden (|1992I ) for MS with m < 30 Mq , repro duce the 



[0/M g] values present in the Galactic Bulge. With the same yields [Gibson 



(119971 ) explains the [0/Fe] values in the intracluster medium and predicts a 
small increase in the C/0 evolu tion at late times for a massive elliptical galaxy. 
By adopting the lArnettI (jl99l[ ) yields in our model we may obtain flatter C/0 
gradients and might not be able to reach the C/0 values o bserved in the H ii 
regions and dwarf stars of the solar vicinity . The yields by lArnettI (|l99l[ ) do 
not consider stellar winds. IChiappini et al.l (|2005l ) using LWY studied the C 
and O evolution in the solar vicinity and the Galactic disk, they reproduce 
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also the C/0 vs 0/H behavior in the solar vicinity, but they predict flatter 
C/0 gradients than the observed ones and they cannot match the high C/0 
values shown by the metal rich star in the solar neighborhood. 

Recently, iMcWilliam et alj (j2007l ) suggest e d tha t the strong metallicity- 
dependent yields for massive stars bv iMaederl (jl992l ). can explain the 0/Mg 
vs 0-Mg/H and the 0/Fe vs Fe/H relations in the Galactic bulge and in 
the solar vicinity, in agreement with our results that favor the HWY model 
over the LWY model. The HWY model includes low O yields at high Z 
and therefore explains: a) the small O increase in the solar vicinity from the 
time the Sun was formed until the present, and b) the flattening of the O 
gradient in the direction of the Galactic center. Massive stars with high Z 
values have strong winds, lose a considerable amount of C and produce high 
C yields. With this C lost, the stars keep a small amount of C needed to 
produce O and consequently their O yields are low. Therefore the C yields 
for massive stars are important at high metallicitics, while their O yields 
are more impo r tant at lo w metallicities, as has been shown previously by 



Akerman et al 



( 2004 ) and Carigi et al. ( ioosl ). By a dding our previous results 



lAkerman et allboO-d : ICarigi et al.ll2005h to those of iMcWilliam et all |2007t ) 



and of this paper we insist that the stellar winds with a high mass loss rate 
are essential to reproduce the high C/O values observed in the disk stars of 
the solar vicinity. 

In the upper panel of Figure 4 we present the evolution of the model that 
assumes HWY and Af^p = 80 M0at different galactocentric radii (4, 8, 12, 
and 16 kpc) that correspond to different final metallicitics [Z ~ 0.028, 0.016, 
0.009, and 0.004, respectively). The AF/AO increase at O > 4 x 10"^ present 
in Figure 4 is due to the lower O yields for massive stars with Z = 0.02. In 
the upper half of Table 1 we show the present-day O values and the AF/ AX^ 
values for each galactocentric radius. We note that ISY j ^Z increases slightly 
with Z for large r values or low O values, while ISY j IS.Z increases significantly 
with Z for short r values or high O values. 

In the lower panel of Figure 4 we show the results for the model with 
LWY and A/„p = 60 M0. This model does not predict an increasing Ay/ AO 
value with increasing O for high O values. In the lower half of Table 1 we 
show the present-day O values and the AF/AX^ values for each galactocentric 
radius. With the LWY model we find higher final O values for lower r values 
because the O yields are higher than than those of the HWY model at high Z. 
Even if the HWY and LWY yields arc identical for low Z, we get lower final O 
values for higher r values for the LWY model because it does not include stars 
with m > 60 Mq. The increase or decrease of O is reflected on the Ay/ AO 
and AY/ /S.Z values because the final Y values are nearly independent of the 
models. 

The helium to oxygen mass ratio, Ay/AO, is an important constraint in 
the study of the chemical evolution of galaxies. We have studied the variation 
of Ay/ AO as O increases in the HWY and LWY evolution models (see Figure 
4). For the HWY model, the model that fits the C/O gradient, we have found 
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Fig. 4. Evolution of Helium vs Oxygen for the Galactic disk at four galactocentric 
distances, 16, 12, 8, and 4 kpc (continuous, long-dashed, short-dashed, and dotted 
lines, respectively). Upper panel: The chemical evolution model assumes Mup ~ 
80 M0and high wind yields. Lower panel: The chemical evolution model assumes 
Muv = 60 M^and low wind yields. See Table 1. Low metallicity H ll regions from 
Peimbert et al.l (120071 1 {filled triangles) and adopting Yp = 0.2477. 



the following relations between Y and O: 



Y ^Yp + AY = Yp + (3.3 ± 0.7)0, 



(1) 



for O < 4.3 X 10~3, and 

Y ^Yp + (3.3 ± 0.7)0 + (0.016 ± 0.003)(O/4.3 x 10"^ - 1)^ (2) 

for 4.3 X 10-3 < O < 9 X lO'^. 

For the LWY model with M„p = 60 MQwe have found the following relation 
between Y and O: 



r = yp + (4.0 ±0.7)0, 



(3) 
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for < O < 11 X 10-3. 



Jimenez et alJ (j2003[ ) from a set of isochrones a nd observations of nearby 



K dwarf stars found that AF/AZ = 2.1 ± 0.4. ICasagrande et al.l (|2007n 
found also that AF/AZ = 2.1 ± 0.9 from the newly computed set of Padoya 
isochrones and obscryations of nearby K dwarf stars. These observational 
results are in yery good agreement with the models presented in Table 1. 

Assuming yield s with low m ass-loss rate due to stellar winds (similar to 
that con sidered by 



Hirschi 



20071) ■ a n d yields without mass loss due to stellar 



winds by Woosley fc Weaver ( 1995 ) , Chiappini et al. ( 20031 ) find for the solar 
vicinity Ay/AZ ^ 2.4 and 1.5, respectively. 

Since stellar winds change significantly the C and O yields, but not the 
Y yields, we are interested to quantify the evolution of the Y contribution 
due to massive stars and due to low and intermediate mass stars at different 
metallicities. For that reason we show in Figure 5 the accumulative percentage 
of Y for four galactocentric distances due to MS and LIMS obtained from the 
HWY model, our successful model for Galactic disk. 

The fraction of helium in the ISM due to MS and and to LIMS depends 
strongly on time, but not on galactocentric radius or Z . At present about half 
of the Ay in the ISM has been produced by MS and half by LIMS. 

The strong dependency on time of the Ay contribution is due to the 
lifetime of the stars. The Ay contribution of LIMS decreases less than 5 
% from 4 to 16 kpc due mainly to the star formation history. In the insidc- 
outside scenario the SFR at 4 kpc is more intense and it is also higher at earlier 
times than at larger distances (see the middle panel in Figure 2), therefore 
the big number of LIMS formed in the first Gyrs at 4 kpc enrich the gas at 
later times. This fact produces the O dilution that can be seen in the right 
hand side panels of Figure 5, just after the model at 4 kpc reaches the value 
of O = 3 X 10^'^ for the first time. 



3. THE HELIUM AND OXYGEN ABUNDANCES OF M17, A HIGH 
METALLICITY GALACTIC H II REGION 

We will compare the predictions of the HWY and the LWY models with 
observations of Y and O in the Galactic disk. At present the best H ii region 
in the Galaxy to derive the Y and O abundances is M17. The reason is that 
the correction for the presence of neutral helium in the abundance determina- 
tion is the smallest for the well obs erved Galactic H ii regions. This is due to 



the high ionization degre e of M17 (jPeimbert et al 



1992 



Esteban et al 



199 



i 



Garcia-Rojas et al.ll2007^ . Due to the large amount of neutral helium present 
in the other well observed Galactic H ii regions, the error in the Y determi- 
nation is at least two times larger than the error for M 17, therefore the Y 
determinations for the other Galactic H ii regions will not be considered in 
this paper. 

To determine very accurate He/H values of a given H ii region we need to 
consider its ionization structure. For objects of low degree of ionization it is 
necessary to consider the presence of He° inside the H"'" zone, while for objects 
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Fig. 5. Cumulative percentage of He as a function of time and oxygen in the ISM, 
produced and ejected by massive stars (MS) and low and intermediate mass stars 
(LIMS) at four galactocentric distances on the Galactic disk (lines as Figure 2). 
The model assumes Mup = 80 M0and high mass loss due to stellar winds, the HWY 
model. 



of high degree of ionization it is necessary to con s ider th e possible presence of a 
He"'""'" zone inside the H+ zone. Peimbert et al.l ( 19921 ). hereafter PTR, found 



for M17 an upper Umit of N(He'^~^)/N{'H~^) of 8 xlO~^ a neghgible amount; 
alternatively they found differences with position of the 7V(He+)/A^H("'") ratio 
correlated with the sulphur ionization structure, result that implies that M17 
is ionization bounded and the presence of a small but non negligible amount 
of Hc° inside the H"'" zone. Therefore for this object the helium abundance is 
given by 

N{Ue) _ N{lle°) + 7V(Hc+) 

iV(H) iV(H+) ■ ^' 

To minimize the effect of the correction for neutral helium we took into 
account only regions M17-1, M17-2, and M17-3, from now on M17 -123, that 



show the highest degree of ionization of the observed regions by PTR. lEsteban et al 
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(|l999f ). and lGarcfa-Roias et alj (|2007| ). as well as the highest accuracy in the 
line intensity determinations by PTR. Following PTR wc will assume that He 
is neutral in the regions where S is once ionized, that is 



iV(Hc") _ N{S+) 
Ar(Hc) ^ N{S) ' 



therefore 
7V(He) 



ICF(Hc) X 



iV(Hc+) 
iV(H+) 



N{S+) 



N{S)'N{S+) 



iV(Hc+) 
7V(H+) ■ 



(5) 



(6) 



To estimate the ICF(He) value we recomputed the N{S'^)/ N{S) ratios derived 
by PTR taking into account that the [S ii] AA 4069 + 4076 lines are blended 
with the O ii A A 4069 and 4076 lines, the correction dimi nishes the [S li] elec- 



tron temperatures from about 12 000 K to about 7700 K (jGarcia-Roias et al 



20071 ). the lower temperatures increase the N (S"'')/iV(S) ratios, and the av- 
erage /CF(He) for M17-123 amounts to 1.035. 

To obtain the N{Hc+)/ N}i{+) value for M17-123 we decided to recom- 
pute the determinations by PTR bas ed on their li n e inte nsities and the new 
helium recombination coeffi cients by Porter et al.l (j2005l ). with the interpo- 



lation formula provided by iPorter et al. (2007). In addition we used the 



hydrogen recombination coefficients by Storev fc Hummer ( 19951) . and the 
collisional contribut i on to the He i lines by Sawev fc Berrington ( 1993h and 
Kingdon fc Ferlandl |l995l ). The optic al depth effects in the t riplet lines were 



estimated from the computations by iBeniamin et al.l (|2002r) . At the tem- 
peratures present in M17 the collisional excitation of the hydrogen lines is 
negligible and was not taken into account. 

To determine the 7V(He"'")/A^(H"'") value we took into account the following 
He I lines AA 3889, 4026, 4471, 4922, 5876, 6678, and 7065. We corrected the 
4922 line intensity by considering that it was blended with the [Fe ill] 4924 
line and t hat the contribution of the Fe line amounted to 5% of the total line 
intensity |Esteban et al.lll999t iGarcia-Roias et al.ll2007l) . The M17-123 line 
intensities adopted are presented in Table 2. 

We did not correct the H and He line intensities for underlying absorp- 
tion, the reasons are the following: a) the average observed equivalent width in 
emission of H(/3), EWe„, (Rl3), amounts to 668 A b) th e predicted EWemiUP) 
for Te = 7000 K amounts to about 2000 A |Allei]ll984h . therefore about 1/3 of 
the continuum is due to the nebular contribution and 2/3 to the dust scattered 
light from OB stars, consequently the underlying stellar absorption only affects 
two thirds of the observed continuum, d) consideri ng the nebular contribu 



tion to the o bservations and based on the models by iGonzalez Delgado et al 



( 1999 ^ 2005 ) for a m o del w ith an age of 2 Myrs as well as the observations by 
Leone fc Lanzafame ( 1998) for A 7065, (since A 7065 was not included in the 
models by Gonzalez Delgado et al. 1999 . 20051 ). we estimated that the EWab 
of the AA 3889, 4026, 4471, 4922, 5876, 6678, and 7065 lines amount to 0.4, 
0.4, 0.4, 0.1, 0.1, 0.2 A respectively, an almost negligible amount considering 
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TABLE 2 

HE I LINE INTENSITIES RELATIVE TO H/3 FOR M17-123 



He I Line 


/ 


3889 


0.1738 ±0.0086 


4026 


0.0240 ±0.0012 


4471 


0.0508 ±0.0013 


4922 


0.0127 ±0.0010 


5876 


0.1549 ±0.0038 


6678 


0.0395 ±0.0010 


7065 


0.0499 ±0.0025 


7281 


0.0064 ±0.0007 



TABLE 3 

PHYSICAL CONDITIONS AND CHEMICAL ABUNDANCES IN M17 



Parameter 


= 0.000 


= 0.036 ±0.013 


T[0 II + III] 


8300 ± 200 


8300 ± 200 


n 


691 ± 246 


744 ± 247 


^3889 


9.5 ±0.8 


11.0 ±0.9 


iV(He+)/iV(H+) 


0.1014 ±0.0014 


0.0982 ±0.0019 


/Ci^(Hc) 


1.035 ±0.010 


1.035 ±0.010 


iV(He)/A^(H) 


0.1049 ±0.0017 


0.1016 ±0.0022 


Y 


0.2926 ±0.0034 


0.2837 ±0.0044 


AF 


0.0403 ± 0.0044^ 


0.0360 ± 0.0053^' 


O 


0.00446 ± 0.00045 


0.00811 ±0.00081 


AF/AO 


9.04 ± 1.35^ 


4.44 ± 0.79 


Z 


0.0101 ±0.0015 


0.0183 ±0.0027 




4.00 ±0.75^ 


1.97 ±0.41 



^Where we adopte d Y„ = 0.2523 ± 0.002 7 
for = 0.000 from lPeimbert et all (|2007h . 

Where we adopte d Y-p = 0.2477 ± 0.002 9 
for / 0.000 from lPeimbert et al.1 (|2007l ). 



the large EWem observed values (see Table 4 in Peimbert et al. 19921 ). e) the 
weighted increase in the helium line intensities amounts to about 0.7%, again 
an almost negligible amount, f) the Balmer lines also show underlying absorp- 
tion and the average correction to the line intensities amoimts to about 0.5%, 
again a negligible amount that cancels to a first approximation the underlying 
correction effect on the He/H line ratios. 
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TABLE 4 

PRESENT DAY VALUES FROM THE GALACTIC DISK MODELS FOR 







R = 


6.75 KPC 




Model 




O(t/w)(10^ 


-3) Ay/ AC AY/AO AY/AFe 


AY/AZ 


HWY A4p 
LWY M„p 


80 Mq 
= 60 Mq 


7.16 
8.28 


6.54 4.47 17.04 
7.26 3.73 17.04 


1.70 
1.62 



To determine the helium physical conditions of the nebula simultaneously 
with the iV(He~'")/j V(H~'") value we used t he maximum likelihood implemen- 
tation presented by Peimbert et al. ( 2002h . This implementation requires as 
inputs: a) the oxygen temperatures, r[0 ill] and r[0 ii], and the oxygen 
ionization degree that provide us with the following restriction 



T[0 II + O III] 



7V(0+)T[0 ii] + N{0++)T[0 III 



(7) 



iV(0+) + A^(0++) 

and b) a large set of helium to hydrogen line intensity ratios. In addition an 
estimate of the electron density, n, in the region where the He lines originate 
is not required but it is useful. This implementation can determine the con- 
ditions of the H II regions either with the restriction of uniform temperature, 
or relaxing this restriction. 

From PTR we adopted r[0 ill] = 8200 ± 20 K. From the 7(3727 /7325) 
ratios for M17-123 by PTR and for M17-3 bv lEsteban et all jl999l ). after 



corre cting the A 7325 A lines for the recombination contribution ( Liu et al 



20001 ). we obtained 8100 ± 1300 K and 9900 ± 1300 K respectively, therefore 
we adopted for r[0 ii] a value of 9000 ± 1000 K. From the T[0 ill] and r[0 ii] 
values and the observations by PTR we find that N{0+)/N(H+) = 0.12 and 
A^(0"''+)/A^(H+) — 0.88. Finally from the previous results and equation (7) 
we obtained that r[0 n + O ill] = 8300 ± 200 K. 

To estimate the electron density we used three determinations the n[S ii] 
for M17-123 by PTR th at amounts to 720 ± 250 cm'^^ the n[0 ii] for M17-3 
bv lEsteban et all (|l999r ) that amounts to 790 ±250 cm^^, and the n[Cl ill] for 
M17-123 that we estimated from the o bserved 1(5518) //(55 38) ratios by PTR 
and the atomic physics parameters bv lKeenan et al.l (|2000l ) that amounts to 
650 ± 450 cm~'^. From the average of these three determinations we adopted 
a value of n = 740 ± 250 cm'^ for M17-123. 

By using as inputs for the maximum likelihood method T[0 ii + O ill] 
= 8300 ± 200 K, n = 740 ± 250 cm^'^, and the helium line intensities pre- 
sented in Table 2 we obtain for M17-123 the n, Tggsg, and 7V(He+)/iV(H+) 
values presented in Table 3. The results for M17-123 are presented in Ta- 
ble 3 for P = 0.000 (constant temperature over the observed volume) and 
for P ^ 0.000 (the temperature variations method). Without the restriction 
of uniform temperature the maximum likelihood of the temperature fluctu- 
ation parameter amounts to P = 0.036 ± 0.013. This P value is in good 
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Fig. 6. Predicted evolution of He vs Oxygen for r = 6.75 kpc. Continuous line: 
model assumes M^p = 80 MQand high wind yield. Dotted line: model assumes 
60 MQand low wind yield. M17 H ll region at r = 6.75 kpc {filled circle, 



0.000) {open circle, = 0.00) 



aKreement with those for M17 derived by PTR, lEsteban et ahl (|l999l ). and 
Garcia-Roias et al.l (|2007l) that are in the 0.033 to 0.045 range; these values 
were determined with two different methods: a) combining the temperature 
derived from the ratio of the Balmer continuum to the Balmer line intensi- 
ties with the temperature derived from /(4363)//(5007) [O ill] ratio, and b) 
combining the O ii recombination line intensities with the A 5007 [O ill] line 
intensities. 

From the mean values of iV(He+) /iV(H+) given in Table 3 and the ICF(He) 
given by equation (6) we obtain Af(He)/7V(H) ratios for M17-123 of 0.1049 
and 0.1016 for = 0.000 and = 0.036 respectively. These values arc similar 
to but more precise than those derived by PTR for M17-123; that amount to 
0.106 and 0.100 for = o.OOO and = 0.040 respectively. 

We obtained the AF/AO and the AF/AZ values presented in Table 3 
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based on the Yp determinations bv lPeimbert et alj (|2007l ). The O abundance 
presented in Table 3 includes both the gaseous and the dus t contribution 



and co rresp onds to the average o f the v alues derived by PTR, lEsteban et al. 



(|l999l ). and iGarcfa-Roias et alJ (j2007l ). these three values are in excellent 
agreement. 

M17 is located at a galactocentric distance of 6.75 kpc, under t he assump 



tion that the Sun is located at a galactocentric distance of 8 kpc (jPias et al 
2002h . In Figure 6 we have plotted the AY/ AO value for t"^ = 0.036, from 



this figure it can be noted that this value is in good agreement, at about the 
one a level, with the Galactic chemical evolution model based on the HWY 
set, and that the O value corresponds to the prediction by the models for a 
galactocentric distance of 6.75 kpc (sec also Table 4); alternatively the M17 
AY/ AO value for = 0.000 is considerably higher, by more than 3cr, than 
the value predicted by the HWY model. Similarly in Figure 6 we compare 
the M17 results with the LWY model, again the values for P = 0.036 are 
in good agreement with the Galactic chemical evolution model while the val- 
ues for P = 0.000 are not. Furthermore from Table 4 it is also found that 
the — 0.036 value for AY/AZ is within one a of the models predictions, 
while the t'^ = 0.000 value for AY/AZ is about 3a away form the models 
predictions. 

With the present accuracy of the AY/ AO and AY/AZ determinations in 
Galactic H ii regions it is not possible to distinguish between the HWY and 
the LWY models. 



4. COMPARISON OF STELLAR AND NEBULAR ABUNDANCES IN 

ORION 

To test the nebular abundances derived with different values for the 
Orion nebula we decided to com pare them with those derived for B star abun- 
dances of the Orion association. ICunha et al. (2006) obtained 12 + log 0/H 



8.70±0.09from 11 B stars, while 



Lanz et al.l (|2008[ ) obtained 12 log Ar/H 



= 6.66 ± 0.06 from 10 B stars. The se results are in excel lent agreement with 
the nebular abundances derived by lEsteban et al.l ([200J) for t"^ ^ 0.00 that 
amount to 8.73 ±0.03 and 6.62 ±0.05 for O and Ar respectively. Alternatively 
the values derived for ^ = 0.00 amount to 8.59 ± 0.03 and 6.50 ± 0.05 for 
O and Ar respectively, values that are about \<j and 3(t smaller than those 
derived from B stars. 



5. ABSOLUTE CALIBRATION OF THE O ABUNDANCE IN THE 

SOLAR VICINITY 

The predicted abundances by chemical evolution models are often com- 
pared with stars and with H ii regions abundances. The most popular com- 
parisons are made with the solar and the Orion nebula abundances. The 
comparisons among the models, the Sun, and the Orion nebula are based on 
the absolute abundances, therefore it is necessary to estimate not only the 
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statistical errors but also the systematic ones in the observational determina- 
tions. 

We will start by considering the solar 0/H abundance. What we want 
from the Sun is the 0/H value when it was formed, the so called initial value, 
and to keep in mind that it is representative of the ISM 4.5 Gyr ago when 
the Sun was formed. We have also to consider that the photospheric and the 
interior solar abundances might be different due to diffusion and gravitational 
settling. 

In Table 5 we present the most popular O/H photospheric determina- 
tions of the last fifty years, the quoted values and the errors are the orig- 
inal ones published by the authors. By looking at the differences among 
the different determinations it is clear that for many determinations prob- 
ably the errors represent only the statistical errors and that systematic er- 
rors have not been taken into account, in sh ort that the to t al err ors have 
been underestimated . The determinations, by Asplund et al. ( 20051 ) and by 



Allende Prietd (|2007l ). are qualitatively different to the previous five because 



they arc based on 3D models while the other five are b ased on ID models 
The last two determinations included in Table 5, those by Caffau et al. (200; 



and Centeno fc Socas-Navarrol ( 2008 ). indicate that the possibility of a further 
revision of the solar abundance is still open. 

The abundances inferred from interior models of the Sun, that are based on 
stereosismological data, are in disagreement with the 3D photospheric models 
and predict heavy elem ent abundances abo ut 0.2 dex higher than the 3D 
photospheric ones (e. g. iBasu fc Antia 1120071 and references therein). 

To compare with our mode ls we will use a s the low O/H value the 3D 
photospheric determination by lAsplund et al. | (|2005r) and as the high O /H 
value the ID photospheric determination bv iGrevesse fc Sauvall (jl998l ). that 
is in good agreement with the helioseismic determination. The next step is 
to have reliable stellar interior models to determine the initia l 0/H value, 
for this purpose we will use the models by Bahcall et al. ( 20061 ) presented in 
Table 6. 

The models by Bahcall et al. ( 20061 ) indicate that the photospheric values 
of Z and Y do not represent the initial values due to diffusion and gravitational 
settling. By assuming that the OjZ ratio is not affected by these processes 
the initial 0/H values correspond to 12 -f log 0/H = 8.70 and to 8.89 for the 
AGS05 and the GS98 surface abundances respectively. To obtain the present 
day ISM value we have to consider that the Sun was formed 4.5 Gyr ago and 
according to the Galactic chemical evolution model by Carigi et al. ( 20051 ) the 
0/H ratio in the ISM has increased by 0.13 dex since the Sun was formed. 
Therefore our determinations of the present day ISM 12 -f 0/H values based 
on the AGS05 and the GS98 abundances amount to 8.83 and 9.02 respectively. 

In Table 5 we also present the most popular 0/H determinations for the 
Orion nebula including the errors presented in the original papers. The pre- 
dictions from the solar abundances and the chemical evolution models have to 
be compared with the total abundances in the nebula that have to include gas 
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TABLE 5 
12 + LOG(0/H) 



Orion nebula ^ 









-/ n nn 

L u.uu 


/2 — n on 




o.yo 






Q 70 4- n 1 






8.77 ± 0.05 


1968(2) 




8.75 ± 0.10 


8.52 ± 0.10 


1977(12) 


8.84 ± 0.07 


1976(3) 






8.49 ± 0.08 


1992(13) 


8.93 ± 0.035 


1989('*) 




8.72 ± 0.07 


8.55 ± 0.07 


1998(14) 


8.83 ± 0.06 


1998(5) 






8.51 ± 0.08 


2000(15) 


8.736± 0.078 


2001(6) 






8.49 ± 0.06 


2003(1^^) 


8.66 ± 0.05 


2005(^) 




8.73 ± 0.03 


8.59 ± 0.03 


2004(") 


8.65 ± 0.03 


2007(8) 










8.86 ± 0.07 


2008(9) 










8.76 ± 0.07 


2008(1°) 










1 Goldbere et al. 


(1960l). 


2 








LambertI (1968). 


3 Ross & AUei 





( 197g) 



Grevesse fc Sauvaj 
(l200lj} 



Holwegei 

hood). 



Anders fc Grevesse 
( 19981). 



( 1989 



Asplund et al 
Allende Prietol JgOOTt ). I 



Centeno fc Socas-Navarro 
CaflFau et all (|2008l ) 

b 



(2008 



10 



11 



Peimbert fc Coster j (Il969l ). 12 



Peimbert fc Torres-Peimbert 



13 



Osterbrock et al.l 



Esteban et all (199 



(| 19771 ). 
( 1992f ). 14 

this value in- 
cludes the fraction of O tied up in dust 
grains that amounts to 0.08 dex. 15 



Deharveng et al.l hood). 16 IPilvugin 



( 20031 ) 



17 



Esteban et all (|2004l ). this 



value includes the fraction of O tied up in 
dust grains that amounts to 0.08 dex. 



and dust. With the exception of the determinations by Esteban et al. ( 19981 



20041 ) that take into account the dust fraction all the other determinations only 



include the gaseous content. The other difference is that there are two possible 
sets of values: a) those that assume constant temperature over the observed 
value given by the 4363/5007 ratio of [ O ill], the = 0.00 case, where is 
the mean square temperature variation l|Peimbertill967l ). or b) those based on 
the O II recombination lines, that are in agreement with those derived from 
the 5007/H/3 ratio taking into account the presence of temperature variations 
over the observed volume, and consequently that P ^ 0.00. 
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TABLE 6 
STANDARD SOLAR MODELS ^ 



Values 


GS98 


AGS05 


initial X 


0.70866 


0.72594 


initial Y 


0.27250 


0.26001 


initial Z 


0.01884 


0.01405 


initial O 


0.00879 


0.00582 


initial 0/H 


8.89 


8.70 


initial AY/AZ 


1.32 


0.88 


initial AY/ AO 


2.82 


2.11 


surface X 


0.7410 


0.7586 


surface Y ^ 


0.2420 ±0.0072 


0.2285 ±0.0067 


surface Z 


0.0170 


0.0122 


surface 0/H 


8.83 ±0.17 


8.66 ±0.17 



^'Stan dard solar models by iBahcall et al.l 
1 20061 ). The GS98 and AGS05 

columns correspond to models with 
the heavy element abundances de- 
rive d from photosph er ic ob servations 



Grevesse fc Sauva] ( 19981 ) and by 



Asplund et al.l (|2005l ) respectively. 

''The error s are the conserva tive ones 
adopted bv lBahcall et al.l (|2006l '). 



Esteban et all (|200J) obtain for the Orion nebula that 12 + 0/H = 8.73 



for ^ 0.00, value that includes the dust correction. By taking into account 
the presenc e of the 0/H Galactic abundance gradient that amounts to —0.044 
dex kpc~^ ( Esteban et al. 2005h we obtain a value of 12 + 0/H = 8.75 for 
the local ISM. This value is smaller than the values estimated for the local 
ISM based on: the solar ph otospheric abundanc es by AGS05 and the GS98, 
the standard solar models bv lBahcall et al.l (j2006[ ). and the chemical evolution 
of the Gal axy that amount to 8 .83 and 9.02 respectively. Similarly from the 
results by lEsteban et al.l (|200J) for the Orion nebula for = 0.00 and the 
observed Galactic gradient we obtain that 12 + 0/H ~ 8.51 for the local 
ISM. From the previous discussion it follows that the best agreement for the 
derived 0/H ISM value is given by the 7^ 0.00 result from Orion and the 
AGSO result from the Sun. From these two determinations we recommend for 
the present day local ISM the value 12 + log 0/H = 8.79 ± 0.08. 

From the previous discussion it follows that the best agreement between 
the solar and the Orion nebula abundances is obtained for the high nebular 
abundances, that are derived from the ^ 0.00 values and that include the 
fraction of atoms tied up in dust grains, and the AGS05 solar value. 
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6. THE SOLAR HELIUM ABUNDANCE 



We want to compare also our Galactic chemical evolution model with the 



helium abundance when the Sun was formed, the initial Y value. Basu fc Antia 



(|2004l ) based on seismic data have derived the Y value in the solar convective 
envelope and amounts to 0.2485 ±0.0034. To derive the initial value we need a 
model of the solar interior that takes into account helium and heavy elements 
diffusion and that agrees with the helium abundance of the en velope. 



Ag ain we have at our disposal the solar interior models bv lBahcall et al 



(|2006l ) presented in Table 6. The GS98 mo del agrees with the Y value in the 
envelope derived bv iBasu fc Antia ( 2004 ). but not with the 0/H value in 
the envelope derived by Asplund et al.l ( 2005 ). On the other h and the AGS05 



model agrees with the O /H value in the enve lope derived by I Asplund et al 
(120051 ) but not with the Y value derived bv iBasu fc Antia I (I2OOJ). Based 
on the discussion of the previous section we conclude that the Orion nebula 
O /H value agrees with the O /H value predic ted by the Ga l actic chemical 
evolution model and the photos pheric value by Asplund et al.l (120051 ) and not 
with the photospheric value by Grevesse fc Sauvall (Il998 ). The discrepancy 
between the photospheric abundances by I Asplund et al.l (|2005l ) and the Y 
value derived from helioseimological data is a very important open problem. 



an exc ellent review discussing this issue has been presented bv lBasu fc Antia 
(I2OO7I). 

Bahcall et al.l (|2006f ) present the initial Y, Z, and O solar values for 
t heir standard solar models, see Table 6. By adopting the Yp value by 



Peimbert et al.l (|2007l ) for ^ o.OO it is also possible to obtain the AY/AZ 



and the AY/ AO values for the GS98 and the AGS05 standard solar models. 
The values so derived are in fair agreement with the predictions of the HWY 
and LWY Galactic chemical evolution models. To try to make a more rigorous 
comparison between the solar interior and Galactic chemical ev olution models 



i t is ne cessary to estimate the errors in the initial solar values bv lBahcall et al 
(j2006f ). a task which is beyond the scope of this paper. 



7. CONCLUSIONS 



Based on the HWY model we find the following equation to estimate the 
initial helium abundance with which stars form in the Galactic disk 



y = Fp + (3.3 ± 0.7)0, 

for O < 4.3 X 10-3, and 

y = Fp + (3.3 ± 0.7)0 + (0.016 ± 0.003)(O/4.3 x 10"^ - l)^ 

for 4.3 X IQ-^ < O < 9 X 10"^ 

The increase of AFe/AZ has to be taken into account in order to deter- 
mine the AY/AZ value based on the [Fe/H] abundances of stars in the solar 
vicinity. 
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High mass loss rates due to stellar winds should be adopted in the evolu- 
tionary stellar models for massive stars of high metallicity, because only the 
Galactic chemical evolution models with HWY can reproduce simultaneously 
the 0/H and C/0 Galactic gradients, the C/O versus O relation in the solar 
vicinity, and the AF/AO value in the inner Galactic disk. 

Based on the O /H value of the Orion nebula and the solar photospheric 
value together with a chemical evolution model of the Galaxy we recommend 
for the present day local ISM a value of 12 + log 0/H = 8.79 ± 0.08, where 
both the gaseous and the dust components of O arc taken into account. 

By comparing the 0/H value of the Orion nebula with the solar value 
we find that the nebular ratio derived using O recombination lines, that is 
equivalent to the use of the forbidden O lines under the adoption of a ^ 0.00, 
is in considerably better agreement with the initial solar value than the Orion 
nebula value derived adopting = 0.00. 



The stellar 0/H an d Ar/H abundance ratios derived by Cunha et al. ( 20061 ) 



and iLanz et al.l (|2008l ) for B stars of the Orion association are in excellent 
agreement with the nebular abundance ratios derived from the ^ 0.00 val- 
ues for the Orion nebula. 

The /^Y/AZ = 1.97 ± 0.41 value derived from observations of M17 for 
~ 0.036 is in very good agreem ent with the 2.1 ±0.4 a n d the 2.1 ±0.9 values 
derived bv I Jimenez et al.l (|2003l ) and ICasagrande et al.l (|2007t ) from K dwarf 
stars of the solar vicinity. On the other hand the value Ay/AZ = 4.00 ± 0.75 
derived from observations of M17 for = 0.000 is not. 

Both Galactic chemical evolution models with the HWY set and the LWY 
set are in agreement with the observed AF/AZ for = 0.036 but not with 
the ISY/ /S.Z for = 0.000. Higher accuracy determinations of AF/AZ for 
high metallicity objects are needed to discriminate between the HWY model 
and the LWY model predictions. 



We arc grateful to Brad Gibson and Antonio Pcimbert for several fruitful 
discussions. We are also grateful to the anonymous referee for some excellent 
suggestions. This work was partly supported by the CONACyT grants 46904 
and 60354. 
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